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Abstract. For each of the eight n-th derivative parameter changing formulas for Gauss 
hypergeometric functions a corresponding fractional integration formula is given. For both 
types of formulas the differential or integral operator is intertwining between two actions 
of the hypergeometric differential operator (for two sets of parameters): a so-called trans¬ 
mutation property. This leads to eight fractional integration formulas and four generalized 
Stieltjes transform formulas for each of the six different explicit solutions of the hypergeo¬ 
metric differential equation, by letting the transforms act on the solutions. By specialization 
two Euler type integral representations for each of the six solutions are obtained. 
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1 Introduction 

This paper has two sources of inspiration. The first aim was to give a complete list of the 
fractional integration formulas corresponding to the eight parameter changing n-th derivative 
formulas for Gauss hypergeometric functions given in [7, 2.8(20)-(27)] and [18, § 15.5]. The 
first fractional generalization of one of these differentiation formulas was given by Bateman [4, 
p. 184]. Fractional generalizations of some further differentiation formulas were given by Askey 
&: Fitch [3, Section 2]. One still missing case was partially handled by Camporesi [5, paragraph 
after (2.28)]. In this paper the full list of the eight fractional integral transformation formulas 
will be given. 

Another observation leading to this paper was that Euler’s integral representation for the 
Gauss hypergeometric function, when written as a fractional integral 

(0 < X < 1, Rec > Re6 > 0), (1-1) 

should have a proof by using the hypergeometric differential equation [1, (2.3.5)] 

“A =0, (1.2) 

*This paper is a contribution to the Special Issue on Exact Solvability and Symmetry Avatars in honour of 
Luc Vinet. The full collection is available at http://www.emis.de/journals/SIGMA/ESSA2014.html 
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where 


La,h,c-,z {f{z)) = {La,b,cf){z) := z{l - z)f{z) + (c - (o + 6 + l)z)f{z) - abf{z), (1.3) 


and that then essentially the same proof should also yield that 


rM 

fix) := - yr\x - dy (1.4) 

J m 


is a solution of the hypergeometric differential equation if m, M and x are as follows: 


(m, M) 

( — 00 , x) 

(-00,0) 

(a;,0) 

(0,x) 

(0,1) 

(0,1) 

(a;, 1) 

(l,a;) 

( 1 , 00 ) 

(x, 00 ) 

X e 

(-00,0) 

(0, 00 ) 

(-00,0) 

(0,1) 

(-00,0) 

( 1 , 00 ) 

(0,1) 

( 1 , 00 ) 

(- 00 ,1) 

( 1 , 00 ) 


This means that m and M in (1.4) are two consecutive points of singularity of the integrand. 
Indeed, we will see that in all listed cases the right-hand side of (1.4) equals a constant multiple 
of one of the six explicit solutions iVj [18, § 15.10(ii)] of the hypergeometric differential equation. 
In fact, as was Sergei Sitnik kindly commenting to me following an earlier version of this paper, 
the above observations (in the fractional integral case where m or M equals x) were already 
made in great detail by Letnikov [14] in 1874 in a paper which seems to have been unobserved 
outside Russia. 

The Euler type integral representations of fractional integral type are specializations of pa¬ 
rameter changing fractional integral transforms acting on some Wj. For instance, (1.1) is a special 
case of Bateman’s fractional integral formula [4, p. 184]: 



ix-y)>" ^ 


dy 


r(c) 


T{c + fi) 


2F1 


a, b 
c + y.' 


(0 < X < 1, Rec > 0, Reu > 0). 


(1.5) 


It will turn out that (1.5), and all other fractional integral formulas for hypergeometric functions 
to be considered, admit a proof by using the hypergeometric differential equation. Here Ta,fe,c) 
given by (1.3), occurs in so-called transmutation formulas, for instance in connection with (1.5): 


d-'a,b,c+fj,',x I X 


d-c-M 


L 




y 


r(h) 


1 (t — 7/^^ ^ 

y^-\La,b,cf)iy) ^ dy 


if G ^^([0,1)), X G (0,1), Rec > 0, Rey> 0). 


( 1 . 6 ) 


Transmutation is a term which occurs in many meanings in science, and even can have various 
meanings in mathematics, but in the sense used here it was first considered in great detail 
by Lions [15], namely as an operator A, often an integral operator, intertwining between two 
differential operators Li and L 2 : 


LiA — AL2- 

In most examples, but not here, Li = (ff jdx^. Then a typical example for L2 would be the 
Bessel type differential operator L2 = df /dx'^ -|- ax~^d/dx. Lions [15] built on earlier work by 
J. Delsarte (1938) and Levitan (1951). Many papers and books on transmutation have appeared 
since then. See the surveys [21, 22] by Sitnik and references given there. 

The cases of (1.4) where m and M are not equal to x are variants of the generalized Stieltjes 
transform, introduced by Widder [26, Section 8], and further developed by many authors, see 
for instance references in [10, 11]. We consider the generalized Stieltjes transform as a transform 
sending f to g oi the form 

l>M 

9ix)= fiy)\y - dy, 

J m 


(1.7) 
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where {m,M) is (—oo,0) or (0,1) or (l,oo), where x is on M outside the integration interval, 
and where the function y i—>■ \y\^~^ f{y) is on the integration interval. The Euler type integral 
representations of this form are special cases of generalized Stieltjes transforms which send a solu¬ 
tion Wi of the hypergeometric differential equation to a solution Wj and change the parameters. 

Some of these formulas can be found in literature, notably in the Bateman project [ 8 ]. There 
[ 8 , 14.4(9)] is essentially the case (m, M) = (—oo,0) of (1.4), while [ 8 , 20.2(10)] is a generalized 
Stieltjes transform sending a 2^1 to a 2 -^ 1 - A formula by Karp & Sitnik [12, Lemma 2] is 
essentially a generalized Stieltjes transform sending a 2^1 to a 3 T 2 , which can be specialized 
to a transform sending 2 E 1 to 2 E 1 . We will give a long list of generalized Stieltjes transforms 
sending some Wi to some Wj, including the two from literature just mentioned. The formulas in 
this list are essentially equivalent: they can all be derived from each other. 

The case ^ = 0 of (1.7) is essentially the classical Stieltjes transform. It will not change the 
parameters of the hypergeometric solutions. A well-known example of this case is [25, (4.61.4)], 
which sends Jacobi polynomials to Jacobi functions of the second kind. 

The idea that formulas for hypergeometric functions can be proved by using the hypergeo¬ 
metric differential equation goes back to Riemann. See Andrews, Askey & Roy [1, Sections 2.3 
and 3.9] how this method can be used for a proof of Pfaff’s and Euler’s transformation formu¬ 
las and of quadratic transformation formulas. Such methods are recently also used by Paris 
and coauthors [13, 24]. They refer to an earlier proof in this spirit by Rainville [19, p. 126] of 
a quadratic transformation formula involving a lEi and a qFi. In this paper the method will be 
applied to integral formulas, but with a different focus. 

Our message is that many formulas for hypergeometric functions have a companion formula 
involving the hypergeometric differential equation, which is more universal because it will imply 
or suggest many formulas involving the various solutions wj of the hypergeometric differential 
equation. A final rigorous proof of these formulas may not use the universal formula (as will be 
often the case in the present paper), but the universal formula is helpful for arriving at these 
formulas and for organizing them. 

Quite probably the ideas of this paper will also work in other situations, for instance for 
Appell hypergeometric series. 

The contents of this paper are as follows. After some preliminaries for hypergeometric func¬ 
tions in Section 2, we illustrate in Section 3 the main ideas of the paper for the special case 
of the Bateman integral (1.5). This takes quite a few pages, but it is less technical than the 
rest of the paper. In Section 4 we state and prove the eight fractional integral transformations 
corresponding to the eight n-th derivative formulas for hypergeometric functions. We give also 
eight corresponding transmutation formulas. In Section 5 we discuss the 48 fractional integration 
formulas for the six solutions Wj, which can be obtained by rewriting the formulas in Section 4. 
Not all formulas will be given explicitly. In Section 6 we give 24 generalized Stieltjes transforms 
sending some Wi to some Wj. In Section 7 we give for each of the six solutions Wj two Euler type 
integral representations. They can all be obtained by specialization of formulas in Sections 5 
and 6 . Einally, Section 8 discusses the connection between generalized Stieltjes transforms of 
different order by fractional integration, and how this leads to connections between formulas in 
Sections 5 and 6 . 


2 Preliminaries about hypergeometric functions 


The Gauss hypergeometric function [7, Chapter 2], [1, Chapter 2], [18, Chapter 15] is defined 
by its power series 


2F1 



F{a,b;c;z) := ^ 
k=0 


{a)k{b)k k 
{c)kk\ 


( 1^1 < 1 ). 


( 2 . 1 ) 
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Here the (complex) parameters a, b, c are taken generically. In particular the series should be 
well-defined (c ^ Z<o); while possibly nicer results in case of terminating series {a or b € Z<o) 
are not paid attention to. The function (2.1) has a one-valued analytic continuation to C\[l, oo). 
It satisfies the hypergeometric differential equation (1.2) and it is uniquely determined as the 
function regular near 0 and equal to 1 at 0 which is annihilated by La^b,c- 

Important transformation formulas, due to Pfaff and Euler, respectively, are [1, (2.2.6), 
(2.2.7)], [18, (15.8.1)]: 

2 F 1 ^ = (1 - z)-“ 2 T 1 ^ (z i [1, 00 )), (2.2) 

2 F 1 ^ = (1 - 2 F 1 ^ [1, 00 )). (2.3) 

Often, a specific formula for F{a,b]C', .) trivially implies another one by the symmetry in a 
and b. For instance, when we will refer to (2.2), we may also mean a similar identity with the 
right-hand side given by (1 — z)~^F{c — a,b]C',zl{z — 1)). An elementary special case of the 
hypergeometric function is 

2Fi(“;_*;2) =(1-2)-“ (2^11.03)). (2.4) 


In the case of generic parameters there are essentially six different explicit solutions of the 
Gauss differential equation Ta,6,c/ = 0 [7, Section 2.9], [18, § 15.10(ii)]. These are one-valued 
analytic functions on the complex plane with a suitable real interval as cut: 


wi{z-, a, b, c) 
W 2 {z; a, b, c) 
W 3 {z;a,b, c) 
W 4 ^{z-, a, b, c) 
W 5 {z-,a,b,c) 
we{z-,a,b,c) 



z^-^ 2 Fi 


a — c-|-l,6 — c-l-1 
2 - c ' 


z-%Fi 

z-\Fi 


/ a,a — c + I 
a — b +1 
/b,b-c+l 
V b — a + 1 ’ 



2 F 1 


a, b 

a + b — c + 1 



(i-zy-'^-^Fi 


c — a,c — b 
c — a — b + 1 



(z^[l,oo)), (2.5) 

(z ^ (- 00 , 0 ] U [ 1 , 00 )), (2.6) 

(z^(-oo,l]), (2.7) 

(z^(-oo,l]), (2.8) 

(z^(-oo,0]), (2.9) 


(2 ^ (— 00 ,0] U [1, 00 )). (2.10) 


Here we had to exclude not just the branch cuts of the 2 Ti’s, but also those of the power factors 
(we assume principal values for the complex powers). In the case of W 2 , W 3 , W 4 , wq we might 
have chosen the branch cuts due to the power factors differently. For instance, a companion 
of W 2 would be 


W 2 {z-, a, b, c) 


{-zf-^2Fr 


f a — c + l,b — 
V 2-c 


c +1 



{z i [0 ,oo)). 


We will also consider the six solutions on subintervals of the real axis including the branch 
cuts of the power factors as below (here we abuse notation by not changing it compared to 
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above) 


wi{x; a, b, c) 
W 2 {x-,a,b, c) 
W 3 {x;a,b, c) 
wa{x; a, b, c) 
w^^x; a, b, c) 
we{x; a, b, c) 





\x\-%Fi 
\x\-\Fi 
2 F 1 


a — c+1,6 — c+1 
2 -c 

a, a — c + 1 
a — 6 + 1 ’ 


\x 


-1 


_i, ^ fb,b- c +I 

- ^ ' ■ X 


b — a + 1 
a, b 


a + b — c + I 


1 — X 


|1 


c — a,c — b 
c — a — b + 1 


; 1 — X 


(xE(-oo,1)), (2.11) 

(xE (-00,0) U (0,1)), (2.12) 

(x E (— 00 , 0) U (1, 00 )), (2.13) 

(x E (— 00 , 0) U (1, 00 )), (2.14) 

(xE(0,(X))), (2.15) 

(xE (0,1) U ( 1 , 00 )). (2.16) 


3 The main ideas illustrated in a special case 


3.1 Transmutation property of a differentiation operator 


In [7, 2.8(20)-(27)] or [18, § 15.5] there is a list of eight parameter changing n-th derivative 
formulas for Gauss hypergeometric functions. One of these is [18, (15.5.4)]: 


2 F 1 


a, b 




1 (c — re) \c — re 


(3.1) 


Formula (3.1) can be proved immediately by power series expansion (2.1). Note also that the 
re-th derivative case follows by iteration of the case re = 1. The case re = 1 can be rewritten as 


Dc-i 2F1 



(c - 1) 2 F 1 


a, b 
c- 1’ 


where 

{Daf)iz) := zf'{z) + af{z) = z"“+^^(z“/( 2 ;)). 


(3.2) 


(3.3) 


Clearly, if / is analytic at 0 then Daf is analytic at 0 and {Daf){0) = o/(0). 

Straightforward computation followed by iteration gives the following transmutation pro¬ 
perty: 


La,b,c—lFc—l — Dc—lLafi^cj (3-^) 

Fa,b,c—nFc—n ' ' ' He—2l^c—1 — n • • • Hc_2Hc_iZ/(j,fe,C' (^■^) 


Formula (3.2) is also a consequence of (3.4). Indeed, by (3.4) La^b,c-i annihilates the left-hand 
side of (3.2). Since this left-hand side is regular at 0, it must be a constant times F{a, 6; c — 1; .) 
with the constant obtained by evaluating both sides at 0. Similarly, (3.1) is a consequence 
of (3.5). 

Because of the above argument, it is natural to consider Dc-iWj{ .‘,a,b, c) not just for j = 1 
but also for the other explicit solutions of the hypergeometric differential equation (j = 2 ,..., 6). 
It will turn out that for all j we get some constant factor times 'Wj{. ;a,b,c — 1). For instance, 
with W 2 given by (2.6), 


(o - c -I- 1)(6 - c -b 1) 


W2{.',a,b,c 


1 ). 


Dc-iW 2 {.;a,b, c) 


2 -c 


(3.6) 
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This is equivalent to the prototypical differentiation formula for the hypergeometric function [4, 
(15.5.4)]: 


dz 


a, b 




_ab 

— — 2.C1 
C 


tt T 1? 5 T 1 
c + 1 




(3.7) 


Corresponding to each of the eight n-th derivative formulas [18, (15.5.2)-(15.5.9)] one can 
write down a transmutation formula like (3.5), and corresponding to each of these transmutation 
formulas one can write down an n-th derivative formula for each of the six solutions Wj. These 
will not be included in the present paper. The transmutation formulas for the first derivative 
operators were earlier given by Derezihski [6, end of Section 3.4]. He calls them commutation 
relations and he relates them to to the Lie algebra so (6) (or equivalently s[(4)). This way 
to associate a Lie algebra (the so-called dynamical symmetry algebra) with a hypergeomet¬ 
ric differential equation was first described by Miller [17]. This was later understood in the 
framework of M-hypergeometric systems, see Saito [20]. 

Remark 3.1. The transformation formulas (2.2) and (2.3) also have a more universal back¬ 
ground because they can be understood from transformation properties of the differential ope¬ 
rator Lafi^c (see also [1, (2.3.10B) and (2.3.10F)]): 

((1 - (^)) = -(1 - , (3.8) 

7.aAc;.((l - = (1 - ^)"-“-'(Lc-a,c-v/)(^). (3.9) 

This suggests identities involving Wi{z;a,b,c), (1 — z)~°'Wj{-^^;a,c — b,c), (1 — z)~^Wk{-^; 
c — a,b, c) and (1 — zY~°‘~^wi{z] c — a,c — b, c). These can indeed be given, but sometimes one 
needs another branch cut for the power factor than chosen in (2.5)-(2.10). 


3.2 Transmutation property of a fractional integral operator 


Both Riemann-Liouville and Weyl fractional integrals will occur in our formulas. See, for 
instance, [8, pp. 181-183] or [16]. 

Formula (3.1) can equivalently be written as a repeated integral, which can be condensed to 
a single integral [1, p. 111]. For this purpose, take 0 < x < 1 and assume Rec > n. We obtain 


2 F 1 



r(c) 

F(c — n) 


1 7? 

y 2x1 



{x — uY ^ 
(n — 1)! 


dy. 


(3.10) 


Note that there are no other terms because, by our assumption, vanishes at y = 0 for 

m = 1,..., n. Formula (3.10) has a fractional extension for n complex with Rec > Ren > 0. In 
rewritten form this is Bateman’s fractional integral formula (1.5), which is often written in the 
form 


1 


t 


c-l 


2 F 1 



ny) 


dt 


frr-M 

T{c + y) Vc + M J 

{z G C\[l, 00), Re c > 0, Re /r > 0). 


(3.11) 


Then (1.5) follows from (3.11) by specializing z to x £ (0,1) and substituting t = y/x. By 
analytic continuation it is sufficient to prove (3.11) for \z\ < 1. There it follows by power series 
expansion (2.1). For (3.11) see also [3, (2.4)], [2, (3.5)] and [1, Theorem 2.2.4, (2.9.6)]. 

Just as we proved (3.2) by (3.4), we can prove (1.5) from the following identity (derived by 
straightforward computation): 

y^~''Li_a,l-b,2-c-,y{y''~^ {x “ yT~^) = La,b,c+ii-,x{x^~''~^{x - 2 /)^"^)- 

We also need (straightforward by integration by parts): 


(3.12) 
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Lemma 3.2. Let f, g G M)). Then, for x G {m, M), 


{La,b,cf){x)g{x) - f{x){Li_a,l-b,2-cg){x) 

= ^(^(1 - x){f'{x)g{x) - f{x)g'{x)) + (c-l + (l-a - b)x) f{x)g{x)). 

Furthermore, let x(l — x)f"{x)g{x), x(l — x)f{x)g”{x), (1 + \x\)f'{x)g{x), (1 + \x\)f{x)g'{x), 
f{x)g{x), as functions ofx, be in L^{{m, M)), and assume that 


- x)f'{x)g{x)\ 

x(l — x)f{x)g'{x) / — ^0 asx 4 ,morxt M. 

(1 + \x\)f{x)g{x) J 

Then 

I'M I'M 

/ {La,b,cf){x)g{x)dx= I f{x){Li_a^i_i,^2-cg){x)dx. (3.13) 

J m J m 

In particular, Li_a^i-b,2-c is the formal adjoint of La^b,c- 
Proof of (1.5) by (3.12). First assume 0 < x < 1, Rec > 1, Re// > 2, We have 
La,b,c+u-,x{x^~''~^ X left-hand side of (1.5)) 

= ^'x,b,c+u-,x{x^~''~^ix - yT~^) dy 

“ Wi {2r(f’f-,y'^'^y'=-Hi-yr-'iy = 0- 

For the first equality sign we twice used a generalized Leibniz rule [9, (1.1)] (allowed because 
the integral in the second line converges absolutely) together with the vanishing of the integrand 
of the left-hand side of (1.5) for y = x (and similarly for the derivative with respect to x of 
that integrand)^. For the second equality apply (3.12) and for the third equality (3.13) (the 
conditions are satisfied). Because 

2 . 1 -C -/1 ^ left-hand side of (1.5) = [ 2F'i( ■,tx'\ — dt, 

Jo \ c J T{fi) 


this is analytic in x at x = 0 , and for x = 0 it takes the value 




r(c) 

F(c-F/x)' 


Thus, by the characterization of the hypergeometric function as solution of the hypergeometric 
differential equation, we have proved (3.12) for Re/i > 2, Rec > 1. These conditions can be 
relaxed by analytic continuation. ■ 


If we replace in the above proof the 2 T 1 by a suitable function / then we obtain the trans¬ 
mutation property ( 1 . 6 ). 

^The involved differentiation under the integral sign got some fame as “Feynman’s trick”, see https://en. 
Wikipedia.org/wiki/Differentiation_under_the_integral_sign#Popular_culture. 
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Remark 3.3. For Re// > 2 (1.6) is equivalent to (3.12). For // = 3, 4, ... (1.6) is also equivalent 
to (3.5). Furthermore, it is sufficient to prove (3.12) for // = 3,4,... because, after division of 
both sides by (x — the two sides depend polynomially on /t. 


Formally, the proof of (1.6) can be extended to showing that 

La,b,c+ti-,x J — dy^ 

It — 7 / 1 ^“^ 

= \yr\La,b,cf){y y dy. (3.14) 

Here m ox M may be equal to x, but not necessarily, and y and x — y should not change sign 
for y G [m,M). This becomes rigorous if /, c and // are such that the first and third equality 
(suitably modified) in the above Proof remain valid in the case of (3.14). In particular, / should 
have suitable vanishing properties at m and M, even stronger if m or M are not equal to 0 or x. 
If (3.14) holds and if Ta,fe,c/ = 0 then 

It — 7 / 1 ^“^ 

La,b,c+t,9 = 0, 5 f(x) := |7/|^~V(l/) — dy. (3.15) 


In this paper many explicit cases of (3.15) will be given, where f = Wi for some i and g = const-tcj 
for some j, with j = i if m or M equals x. As an example it can be derived that 


/ X 

{-yy-^W 2 {y;a,b, c) 

-OO 


(x 


- 2 /)^-^ ^ _ F(a-c-// + l)r(6-c-// + l) r(2 
Yig) ' ^ “ r^(a - c + 1) r(6 - c +1)^ r(2 - 

X {—xy^^~^W2{x; a,b,c + y) 


-c) 
c- y) 
(3.16) 


(x < 0, Re (a — c — /i + 1), Re (6 — c — // + 1), Re /i > 0). 


This is the fractional generalization of the iteration of (3.6). It can be equivalently written as 


f 


2 F 1 


a, b 


(x-y)M-i T{a-y)rib-y) F(c) ^ 
dy = —-^^ 77 ^— 777 - 7 2 x 1 


r(/^) 


F(a) r(6) T{c-y) 

(x <0, Rea, Reft > Re/i > 0). 


a — y,b — y 
c — y ’ 


(3.17) 


This is the fractional generalization of the iteration of (3.7). A different proof will be given in 
the next section. 

Another explicit case of (3.15) which we will meet is 



y^ ^WG{y;a,b,c){y 


xY ^ dy 


F(a — c — y + l)r(ft — c — y + l)F(c — a — 
r{2 - c - y)T{l - y) 

(0 < X < 1, Re (a — c — /i + 1), Re (ft — c — 


i^ 2 (x; a, ft, c +/i) 
y + 1), Re (c — a — ft + 1) > 0). 


(3.18) 


The left-hand side is no longer of fractional integral type, but it is a generalized Stieltjes trans¬ 
form, to which we will return in a moment. 


3.3 Euler type integral representations 


When we replace c, /i by ft, c — ft in (1.5) or (3.11) and use (2.4) then we obtain Euler’s integral 
representation as fractional integral ( 1 . 1 ), or in its most used form 

u 

6—1/1 j.\c—6—1/1 a 


2 El 


a, ft 


5 ; = 


r(c) 




r(ft)r(c-ft) Jo 

{z G C\[l,cx)), Rec > Reft > 0). 


(3.19) 
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The right-hand side of (1.1) is annihilated by La^b,c\x- This is a consequence of the transmutation 
property (1.6), by which La^b,c;x acting on the right-hand side of (1.1) is equal to 


const •X 


l-c 


y’’ ^La,b,b-,y{{'i- - y) ‘")ix-yy ^ ^dy = 0 


because (1 — y) “is annihilated by La^^-y. This last fact is also apparent from 

The more general transmutation property (3.14), considered with c, y, replaced by 6, c — 6, 
suggests that each wj has an Euler type integral representation (1.4). This is indeed the case, 
as we already briefly indicated after (1.4). 

We can write (3.19) also equivalently as a generalized Stieltjes transform 


W3{x;a,b,c) = 


T{a-b + l) 




r(a - c + l)r(c - b) Jq 
(x > 1, Re (a -|- 1) > Re c > Re 6). 


(3.20) 


It will turn out that, more generally and analogous to (1.4), 

rM 

9{x)= / \yr^\l-yr^-^\x-yrdy (3.21) 

J m 

is a solution of the hypergeometric differential equation if m, M and x are as listed after (1.4). 

Although most Euler type integral representations of fractional integral type are equivalent 
to some integral representation of generalized Stieltjes transform type by a change of integration 
variable, this is no longer true for the integral transforms mapping Wi to Wj which specialize to 
an Euler type integral representation. 


3.4 Generalized Stieltjes transforms as transmntation operators 


We already observed in Section 1 that variants of fractional integral transforms and the Euler 
integral representation for hypergeometric functions naturally lead to formulas involving a ge¬ 
neralized Stieltjes transform. In the definition by Widder [26, Section 8] the generalized Stieltjes 
transform sends a suitable measure a or function (j) (with da(t) = (p{t) dt) to a function / analytic 
on C\(—oo, 0]: 


da{t) 

0 {z + t)P 


/(^)- 


The special case p = 1 gives the classical Stieltjes transform. In order to have analytic expressions 
similar to the ones in fractional integral transforms, we will work with transforms (1.7). 

Transforms of generalized Stieltjes type have transmutation properties. Eor instance, from 
(3.14) we have, associated with (3.18), the intertwining property 


L 


a,b,c+tJ. 


■,x(^x^ V(y)(2/-a:)^ ^dy^ 

D 

y''~^iLa,b,cf)iy){y - dy 




(0 < X < 1). 


Noteworthy is the case y = 0. Then we have the same hypergeometric differential operator on 
both sides. The Stieltjes transform will then map solutions of the differential equation to other 
solutions. 
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Karp & Sitnik [12, Lemma 2] proved the following formula: 

(.Ve-V-V ' - 

r(5)r(c)r((i + e — b — c) f a,b,c 


mm 


3 F 2 


d, e ’ 


{z G C\[l, cx)), Re 6, Rec, Re (d + e — 6 — c) >0). 


(3.22) 


If we replace z by z~^ then we recognize the formula as a generalized Stieltjes transform sending 
a 2 F 1 to a 3 L 2 of general parameters. (For more general formulas expressing pFp-i functions 
as generalized Stieltjes transforms see [11] and references given there.) If moreover d = a the 
transform in [19] sends a 2 F 1 to a 2 Fi'. 




a+e—b—c—1 


2F1 


a — c,e — c 
a + e — b — c’ 


; 1 — t) (2 — t) dt 


r(5)r(c)r(a + e- 6 -c) fb,c . 

r(a)r(e) ^ \ e ’ 

{z G C\(— 00 , 1], Re6, Rec, Re (d + e — 6 — c) > 0). 


(3.23) 


All generalized Stieltjes transforms mapping a solution wt to a solution wj we know can be 
obtained from (3.23) by change of parameters, change of integration variable, and application 
of (2.2) and (2.3). 


4 The eight fractional integral transformations 
of the Gauss hypergeometric function 

Each of the eight n-th derivative formulas in [18, § 15.5] has a fractional generalization. Some of 
these are very well-known, but others were hardly known until now. They fall apart into three 
families. Within a family the formulas follow from each other by application of (2.2) or (2.3). 
We will use shorthand names for the eight cases of which the meaning will be obvious. The 
division of the cases into families and their correspondence with the n-th derivative formulas is 
as follows: 


family 

I 

I 

I 

II 

II 

III 

III 

III 

case 

c+ 

0-|-, C-h 

0“l“, b-h, c-{- 

a— 

a+ 

a-, b-, 

c— a—, c- 

c— 

[18, § 15.5] 

(4) 

(8) 

(9) 

(3) 

(5) 

(1) 

(7) 

(6) 


4.1 Family I 

Case c+. We already discussed Bateman’s fractional integral formula (1.5). From (3.11) we 
get another variant of (1.5): 


i-ym 2F1 


a, b 


;y 


{y-xY ^ 

m 


dy = 

r(c-i-/i) \c + yi 

(x < 0 , Rec > 0, Re// > 0), 


\x 


which we can write together with (1.5) in a unified way as: 




' 0<y/x<l 


a, b 




r(//) 


TY + n) 


2 F 1 


a, b 

c + Y 


{x G (—00,0) U (0,1), Rec > 0, Re/i > 0). 


(4.1) 
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Case a+, c+. By using (2.2) in (4.1) we arrive at 
Jo<y/x<l \ ^ ) ^(/^ 


dy 


r(c) 


X 


V C + /i J 


r(c + /i) 

(x G (—oo, 0) U (0,1), Rec > 0, Re// > 0). 


(4.2) 


Case a+, 6+, c+. By using (2.3) in (4.1) we arrive at 

ny) 

'a + //, 6 + /i 


Jo<y/x<l \ ^ J ^(l*) 


dy 


r(c) 


\x 


C+fl—l 


(1 -x) 


a+6—c+/i 


2 F 1 


r(c + //) 

(x G (—00,0) U (0,1), Rec > 0, Re/i > 0). 
See also Askey Sz Fitch [3, (2.11)]. 

4.2 Family II 

Case a—. Askey & Fitch [3, (2.10)] give 


c + /i 


;x 


(4.3) 


t 


a—fi—1 


2 F 1 


c J F(//) F(a) 


a — fi,b 




(z ^ (1,00), Rea > Re/i > 0). 


(4.4) 


The proof is as for (3.11). For \z\ < 1 the formula follows by power series expansion (2.1), 
and next the general case follows by analytic continuation. Formula (4.4) implies the following 
fractional integral formula; 

'a,b \ Jx — y\^~^ 


'0<y/x<l 


|y| 


a—fi—1 


2 F 1 




dy = 


F(a-/i) fa- fi,b 


F(a) 


Fi 


(x G (—00,0) U (0,1), Re a > Re /t > 0). 


(4.5) 


Case a+. By using (2.2) or (2.3) in (4.5) we arrive at 


'o<y/x<i 


|yr-“-'^-i(i-2/)“+'-SFi 


a,b \\x-y\>^^ 


F(c — a — /i) 
F(c — a) 


X 


c—a—1 


(1 -x) 


a+6—c+/i 


2F1 


a + ii,b 


;x 


(x G (—00, 0) U (0,1), Re (c — a) > Re// > 0). 


(4.6) 


4.3 Family III 


Case a—, h—, c—. The following generalizes a formula of Camporesi [5, paragraph af¬ 
ter (2.28)]. 


Proposition 4.1. IFe have 


f 


a,b\{x-yY^^^ T{a-y)T{b-y) F(c) 

y --“ 7 / = —-,^ 77 ^- 777 - 7 2^1 


c ’V F(/i) F(a) F(6) F(c -//) 

(x < 1, Re a, Re 6 > Re /i >0). 


a — /t, 6 — /i 
c — // ’ 


(4.7) 
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Proof. Formula (4.5) can be rewritten as: 




\y\ ’ ■,y 


a,b^ _A \x - y\ 


— 1 




, dy = — ' ' \x 

r(/i) r(a) ' 


(x G (—oo, 0) U (1, oo), Re a > Re/i > 0). 


(4.8) 


Combination with [7, 2.10(2) and (5)] (or with [18, (15.8.2)]) yields (4.7) for x < 0. This, in its 
turn, implies by transformation of integration variable and by analytic continuation that 


2 F 1 


“’^•1 (t - 1)^ ^ _ r(a-y)r(b-/j) r(c) 

c ’ J r(/i) r(a) r( 5 ) r(c-y) 


X (-z)-f^2Fi 


a — y,b — y 
c — y ’ 


1 + z 


(z ^ [0,oo), Rea > Re^ >0). 


(4.9) 


From that we see that (4.7) holds for x < 1. 


Case a—, c—. 

Ai-y) 

J X 

x(l 


By using (2.2) in (4.7) we arrive at 
a_^_i j:. f a,b \ {y-x^-^ r{a-y)r{c-b-y) F(c) 

— a:)““^ 2 r'i| ;a:| (a: < 1, Rea, Re (c — fc) > Re u > 0). 

\ c — u / 


(4.10) 


Case c— . By using (2.3) in (4.7) we arrive at 


r il-yr^^-^2Fi 

J — OO 



(x — yY ^ 

ny) 


dy 


T{c — a — y) T{c — b — y) r(c) 
r(c —a) F(c —6) F(c —^) 


x(l 


xY+b-c+t, 


a, b 
c — y 



(x < 1, Re (c — a). Re (c — b) > Rey > 0). (4.11) 


4.4 Transmutation formulas 

Corresponding to case c+ above, we gave already (3.12), which gives rise to the transmutation 
formula (1.6), and by which a proof of (4.1) can be given. Analogues of (3.12) and (1.6) can be 
given for all cases above. These have the general form 


F 


a\b\c';x 




and 


La',b',c'-,x(^J^f{y)^^^\x-yY ^dy'^ = J^iLa,b,cf)iy) 


v{x)V2{x) 

w{y)w 2 {y) 
v{x)V2{x) 


(4.12) 


x-yY ^dy. (4.13) 


The data to be used in these two formulas for the eight cases are specified in the table below. 
In general, w{y) is a product of powers of \y\ and 1 — y, and W 2 {y) is equal to \y\, 1 — y or 1. 
Similarly for v{x) and U 2 (x), respectively. The variable y in (4.12) ranges over an open interval I 
which is also the integration interval in (4.13). The interval I has endpoint x at one side and 
endpoint xq = 0,1 or —oo at the other side. The variable x in (4.12) and (4.13) ranges over 
some open subset J of M. The function / in (4.13) should be in C'^{J) and should moreover 
satisfy certain growth conditions at xq, to be specified in a moment. The parameter y can be 
arbitrarily complex in (4.12) but should satisfy Re/r > 2 in (4.13). 
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case 

a! 

b' 

c' 

Xo 

w{y)/v{x) 

W2{y)/V2{x) 

J 

c+ 

a 

b 

c -1- y 

0 


1 

(-oo,0)U(0,l) 

a-t, c+ 

a + y 

b 

c -1- y 

0 


1-y 

(-oo,0)U(0,l) 

|x|^+'"-i(l - xf-<^ 

|y|c-l(l_y)a+6-c 

1 — X 

1 

aH“, cH- 

a + y 

b + y 

c -1- y 

0 

(-oo,0)U(0,l) 


a— 

ci-\- 

a — jl 

a-\- fl 

b 

b 

c 

0 

0 

|y|a-M-l/|2.|a-l 
|y|c-a-/i-1^2^ _ yY+b-c 

\y\/\x\ 

\y\/\x\ 

(-oo,0)U(0,l) 

(-oo,0)U(0,l) 


|^|c-a-l(l _ Y^+b-c+y. 

a—, b—, c— 

a — fl 

b — y 

c — fl 

— OO 

1 

1 

(-00,1) 

a—, c— 

a — fl 

b 

c — fl 

1 

(1 — x)°‘~^ 

i-y 

1 — X 

(-00,1) 


a 

h 

c — fl 

— OO 

(1 - = 

1 

(-00,1) 



^2^ _ j.^a+b—C+/J 


Verification of (4.12) in the eight cases is by straightforward computation, possibly using com¬ 
puter algebra. 

In order to find the growth conditions in (4.13) for / at xq we recall that (4.13) is obtained 
from the string of equalities 


L, 


a' ^c'\x 




- f fl \T ( ^(2/)^2(y)| ,a-l\ 


v{x)v2{x) 
w{y)w2{y). 


dy 


Jl v{x)V2{x) 


In the various steps we have to consider the singularities for y at x and at xq. The first singularity 
is already taken into account by the condition Rey >2. As for the singularity at xq the first 
and the third equality need extra assumptions (which will also imply that the four parts of the 
string are well defined). For the first equality we need 


f{y)w{y){l + |y|)^®^ ^ as a function of y is at y = xq. 


For the third equality the needed assumptions follow from Lemma 3.2: 


/"(y)y(l - y)u:(y)u; 2 (y)(l + ^ 

f{y)wiy)w2{y){i + 

/(y)y(l - y)^{w{y)w 2 {y){l + |y|)^®^“^) 

f{y)-^{wiy)w2iy){i + |y|)^®^) 

f{y)w{y)w2{y){l + |y|)^^®^-i 


> as functions of y are at y = xq, 


and 


/'(y)y(l-y)M;(y)ri; 2 (y)(l + |y|)^®^ ^ 
/(y)y(l - y)^{w{y)w 2 {y){l + |y|)^®^"^) 
f{y)w{y)w2{y)il + |y|)^®^ 


0 as y 


XQ. 
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Thus, for given /, the identity (4.13) is settled under certain constraints for a, b, c and //. 
But then the identity will be valid under relaxed constraints on a, b, c and ^ such that both 
sides of (4.13) are analytic in these four parameters, i.e., if Re^u > 2 and 

/"(y)2/(l - y)w{y)w2{y){l + j 

f'{y)w{y)w 2 {y){l + / as functions of y are at y = xq. 

f{y)wiy){'^ + \y\)^''^~^ J 


5 Fractional integral transformations for the six solutions 
of the hypergeometric differential equation 


Formula (4.1) is a fractional integral transformation of type c+ for the solution wi of the 
hypergeometric differential equation. It turns out that for all six solutions wt there is such 
a transformation of the form (. .;a,b,c) —>■ (. .; a, 6 , c + /r). These can all be 

obtained by rewriting fractional integral transformations for 2 ^ 1 ( 0 , b]C] .) of various types given 
in Section 4, namely types c+; a—, b—, c—; a—; a—; c—; a+, 6 +, c+, respectively. Here we list 
these formulas. Each formula is preceded by the formula number of the formula in Section 4 to 
which it reduces: 


(4.1): 


f I ic-i f u ^{x-yY ^ 

/ |y| Wi{y,a,b,c) —z—— 

/0<y/x<l ^ 


dy = 


r(c) 


^wi{x\a,biC + y) 


T{c + y) 

{x G (— 00 , 0) U (0,1), Re c > 0, Re /U > 0); (5-1) 

,. Y I ic-i / T ~,{x-y)^~^ j T{a-c-y + l) 

4.7): / y ^W 2 {y;a,b,c) ^ — dy = — - 

J-00 r(/i) r(a-c + i) 

T{b — c — y + l) r(2 — c) \xY+^^-'^u;y(x■ a b c + n) 

T{b-c + l) T{2-c-yy' W 2 [x,a,b,c + y) 

{x G (— 00 ,1), Re (a — c + 1), Re (6 — c + 1) > Re /r > 0); (5-2) 

(4.8): [ \yY~^W3{y,a,b,c) ^^ — dy = ^^ \xY^'^~^W 3 {x] a, fe, c + y) 

Jy/x>i r(/i) r(a-c+l) 

{x G (— 00 , 0) U (1, 00 ), Re (a — c + 1) > Re/r > 0); (5-3) 

(4.8): f \yY~^W4{y,a,b,c) ^^ — dy = ^ \xY^f^~^W4ix; a,b,c + y) 

Jy/x>i r(/i) r(6-c + l) 

{x G (— 00 , 0) U (1, 00 ), Re (6 — c + 1) > Re /U > 0); (5-4:) 


(4.11): / y^ ^W5{y;a,b,c) 


{y - xY ^ T{b-c-y + l) 
■dy = 


r(;u) r(6-c + l) 

T{a-c-y + l) r(a + 6 -c+l) c+u-i / i. , ^ 

' \ A -x;a,6,c + // 

r(a —c + 1) T[a + b — c — y + l) 

{x G (0, 00 ), Re (6 — c + 1), Re (a — c + 1) > Re (y) > 0); 


(5.5) 


(4.3): 


y" ^wYy;aY,c^ 




\y-x 


iq<Y^<i ' " r (//) 

l — x 


■dy 


^ ^WGix-,a,b,c + y) 

r{^c — d — b -\- fjj ij 

(x G (0, 1) U (1, 00 ), Re (c — a — 6 + 1) > 0, Re ^ > 0). 


(5.6) 


It is a straightforward exercise to list also the fractional integral transformation formulas 
of the Wi corresponding to the other seven types. We do not list all these formulas here, but 
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only give their essential behaviour in the following table. Here an entry in ith row, jth column 
(place (i, j)) tells us that the transformation formula for Wj of type given at (z, 1) can be reduced 
to the transformation formula for wi of type given at (z, j). 


Wl 

W2 

m 

Wi 

W5 

We 

c+ 

a—, b—, c— 

a— 

a— 

c— 

o-f-, b-\-, c-\- 

a.-|-, c-\- 

a—, c— 

0 “l“, c-\- 

a—, c— 

a+ 

a+ 

a.-|-, b-\-, c-\- 

c— 

(2+ 

a+ 

0.-1-, b-\-, c-\- 

c— 

a— 

a— 

a—, b—, c— 

c+ 

a—, c— 

o-j-, c-\- 

a+ 

a+ 

cz-|-, b-\-, c-\- 

c— 

c-\- 

a—, c— 

a—, 6—, c— 

c+ 

a— 

a— 

a—, b—, c— 

c+ 

a—, c— 

0 “l“, c-\- 

a—, c— 

CL-\-, C-\- 

a— 

a— 

c— 

b-\-, c-\- 

a+ 

a+ 

c+ 

a—, b—, c— 


The cases c+; c—; a—, 6—, c—; a+, 6+, c+ occur 5 times, the cases a+, c+; a—, c— 6 times, and 
the cases a—; a+ 8 times. That these numbers are not all equal will be caused by the symmetry 
in a and b of the hypergeometric function. 

Some more examples from the above table which we need later (because of specialization to 
Euler type integral representations) are: 


/ |yr ^W 2 {y,a,b, — dy = ^W 2 ix-,a - y,b,c) 

Jo<v<i r(/z) r(a-c + l) 

X 

(x G (—oo, 0) U (0, 1), Re (a — c + 1) > Re/z > 0), 
f ( u j r(6-/z) 

/ W 4 [y, a, b, c) --- dy = W 4 [x-, a - y,b- y,c- y) 

J3^>l r(/z) r(6) 

X 

{x G (—oo, 0) U (1, oo). Re 6 > Re^ > 0), 

f r I, J T{c-a-b+l) , ^ , 

/ we{y, a, b, c) dy = —-— , , —-zn6(x; a - y,b- y,c- y) 

Jo<^<i r(^) r(c-a-6 + ^ + l) 

(x G (0, 1) U (1, oo). Re (c — a — 6 + 1) > 0, Re /z > 0). 


(5.7) 


(5.8) 


(5.9) 


6 Generalized Stieltjes transforms between solutions 
of the hypergeometric differential equation 


First we recall the proof of (3.22) as given by Karp &: Sitnik [12, Lemma 2]. Expand for | 2 :| < 1 
the left-hand side of (3.22) as 


k=0 


Z- / 2F1 


By the limit case x f 1 of (1.5) this equals 

kl 


{F)k k + e — b — c)r(6 + k) 
2 ^—^ -^^7^--2^1 


k=0 


T{d + e — c + k) 


d — c,e — c 
d + e — b — c’ 


d — c, e — c 
d + e — c + k' 


; t ) dt. 


By the Gauss summation formula [1, Theorem 2.2.2] we get 


E 

k=0 


{a)k fcr((i + e — b — c)r(6 + k) r(c + k)r{d + e — c + k) 


k\ 


r(d + e — c-\- k) 


T{d + k)T{e + k) 


which is equal to the right-hand side of (3.22). All steps can be rigorously justified because of 
the constraints in (3.22). 
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We now list the generalized Stieltjes transforms mapping a solution Wi to a solution Wj. They 
can all be obtained from the special case (3.23) of (3.22) by change of parameters, change of 
integration variable, and application of (2.2) and (2.3). In particular, (6.10) and (6.12) below 
are simple rewritings of (3.23) by a change of parameters. Also, (6.1) below is a rewriting of [8, 
20 . 2 ( 10 )]. 

The formulas in the list are grouped in cases similar to the cases in Section 4. In all formulas 
the constraints are that x is in an interval (—oo, 0), (0,1) or (1, oo) which does not coincide with 
the integration interval, and that the arguments of the three gamma functions in the numerator 
on the right-hand side have positive real parts. 


Case c+. 

rO 


{-yf ^wi{y]a,h,c){x-yY ^ dy 

= ^+m). (^.d 

r(a + b — c — y + l)r(l — n) 

) 

y''~^W6{y; a, b, c){y - xY~^ dy 

Via — c — a + l)T(h — c — u + l)T{c — a — h + 1) , , x / % 

= ^ a,b,c + /i . 6.2 

r {2 - c - n)T[l - fi) 


Case a+, c+. 


/ 2/'" ^Y-yf '^Wi{y;a,b,c)\x - yY ^dy 

Jo 

r(—a — /i + i)r (6 — c — y, + i)r(c) 
r (6 — a — y + l)r(l — y) 

X \x\^^^~^\x — l|^“'^rc4(x; a + y^h,c + y), 

/ OO 

y''~^iy - a, b, c){y - xY~^ dy 

r(—a — y + l)r (6 — c — y + l)r(a -6-1-1) 
T{-c- y + 2 )T{l- y) 

X — l\^~'^W2{x] a + y,h,c + y). 


(6.3) 


(6.4) 


Case 6+, c+. 


/ 2/'" ^Y-yT ^wi{y]a,b,c)\x-yY ^ dy 

Jo 

r(—6 — y + l)r(a — c — y + l)r(c) 
r(a — b — y+ l)r(l — y) 

X — l|““% 3 (x; a,b + yiC + y), 

/ OO 

y''~Yy - a, 6 , c)|x - yY~^ dy 

r(—6 — y + l)r(a — c — y + l)r(6 — a -|-1) 
T{-c- y + 2)TY - y) 

X — l|““%2(a;; a,b + y,c + y). 


(6.5) 


( 6 . 6 ) 
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Case a+, 6+, c+. 

[ a, b, c){x - dy 

J —oo 

_ r(-a - /i + i)r(-6 - /i + i)r(c) 
r(c — a — b — y + l)r(l — y) 

X a + y,b + y,c + y), (6.7) 

/ OO 

y''~^{y - a, b, c){y - xY~^ dy 

r(—a — y + l)r(—6 — y + l)r(a + 6 — c + 1) 
r{2 - c - y)T{l - y) 

X - x)“+^"‘'+'"u; 2 (a:; a + y,b + y, c + y). (6.8) 


Case a- 


j (-y)“ ^ ^wYy,a,b,c){x -yY ^ dy 

J —OO 




/ 2/“ ^ ^W6{y;a,b,c)\x - y\>^ ^ dy 
Jo 


r{a-c-y + l)r(a - y)r{c - a - b + 1 ) , .g-i . . ^ ^ 

r(a-6-^ + l)r(l-;u) ' ' 


(6.9) 


( 6 . 10 ) 


Case b- 


/ {-yY ^ ^W3{y]a,b,c){x-yY ^ dy 

J — OO 

= r(6 - y)nb -c-y + l)r(a - b + Y b-i . « 5 _ ) 

r{a + b-c-y + l)r{l-y) w^{x,a,b y,c), 

j y^-^-^wYy,a,b,c)\x-yY-Uy 

= r(fe - c - + l)r(b - y)Tic -a-b+1 ) , , 6-1 . b - u c) 

T{b-a-y + l)T{l-y) ' ' wYx,aY y,c). 


( 6 . 11 ) 


( 6 . 12 ) 


Case a+. 


i-y) 


c—a—fi—1 


(l-y)“+'' ''w3{y,a,b,c){x - yY ^ dy 


(6.13) 


r(-a -y + l)r(c - g - ^)r(a - b + 1) c-a-u^ _ .a+b-c+, . ^ ^ . 

T{c-a-h-y + l)T{l-y) w^[^x,a + y,yc), 


y^-^-^-\l - yT^'^-^w^iy a, 6, c)\x - yr^ dy 
Jo 

r(—a — y + l)r(c — a — y)T{a + 6 — c + 1) 
r(6 — a — y + l)r(l — y) 

X |a:|^“““^|l — xY~^^~‘^~^W 4 {x; a + y,b,c). 


(6.14) 
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Case b+ 

fO 


[ - y)“+'-%4(y;«, b, c)(x - dy 

J — oo 

T{—b — /U + l)r(c — b — /u)r (6 — a + 1 ) 
r(c — a — b — fi + l)r(l — y,) 

X — xy^'^~^~^^WQ{x; a,b + fi, c), 


c-b-ii-l(-i _ „,\a+b-c. 


{1 - y)-+'’-W5iy,a,b,c)\x - dy 
T{—b — /U + l)r(c — b — /u)r(a + 6 — c + 1) 

r(a — b — fi + i)r(i — 

X — x\°'~^’^~'^~^W3{x] a,b + fi, c). 


Case a—, b—, c—. 

/■° 

/ W 2 {y,a,b,c){x-yY ^ dy 

J —OO 

T{a-n)T{b-fi)T{2-c) 
r(a + b — c — y, + l)r(l — y) 

/ OO 

we{y;a,b,c){y - xY~^ dy 

r(a — y)T{b — y)T{c — a — b + 1) 
T{c - y)r{l - y) 


W 5 {x; a — y,b — y,c — y), 


wi{x] a — y, b — y, c — y). 


(6.15) 


(6.16) 


(6.17) 


(6.18) 


Case a—, c- 


/ (l-y)“ ^ ^W 2 {y,a,b,c)\x-y\^^ ^ dy 
Jo 


T{c-b-y)r{a-y)r{2-c) _i 

'x-l\ ws{x;a- y,b,c- y), 


r(a — b — y + l)r(l — y) 

/ OO 

{y - l a, b, c){y - xY~^ dy 


T{c- y)r{l - y) 


(6.19) 


( 6 . 20 ) 


Case b—, c—. 


/ (1-y)^ ^ ^W 2 iy,a,b,c)\x -y\>^ ^ dy 

Jo 

T{c-a-y)T{b-y)T{2-c) 

/ OO 

{y - l)^“^“^u> 3 (y; a, b, c){y - xY~^ dy 
r(c — a — y)T{b — y)T{a — 6+1) 


( 6 . 21 ) 


r(c-//)r(i -^) 


\1 — xy ^wi{x;a,b — y,c — y). 


( 6 . 22 ) 
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Case c- 


[ (1-2/)“'^'’ ‘'w 2 {y]a,b,c){x-vY ^ dy 

J — OO 


r(c — a — b — y. + l)r(l — y) 


(6.23) 


/ OO 

{y - a, b, c){y - xY~^ dy 


r(c - a - y)r{c -b- y)r{a + b-c+ 1) 


T{c- y)T{l- y) 


'Wi{x] a,b, c — y). (6.24) 


We summarize the results of the above list in the following table. In the box in row Wi and 
column Wj the type is given of the generalized Stieltjes transform sending Wi to Wj. 



Wl 

W2 

W3 

Wi 

W3 

We 

Wl 



6 +, c+ 

a+, c+ 

C+ 

a+, 6+, c+ 

W2 



a—, c— 

6 —, c— 

a—, b—, c— 

c— 

W3 

b—, c— 

fl+, c+ 



b- 


W4 

a—, c— 

5T, c+ 



a— 

b+ 

W5 

c— 

o+, 6+, c+ 

b+ 

a+ 



We 

a—, b—, c— 

c+ 

a— 

b- 




Due to the denominator factor r(l — /x) all generalized Stieltjes transforms above become zero 
if ;U is a positive integer satisfying the constraints. 


7 Euler type integral representations 

Below we give the explicit Euler type integral representations which have form (1.4) or (3.21). 
Euler type integrals ^vith integrand |1 — y\~°'\x — 


wYx;a,b,c) = [ \y\^-Hl-y)--\x-yr^-Uy 

T{b)T[c-b) Jo<y/x<i 

{x G (—OO, 0) U (0, 1), Re c > Re 6 > 0), (7.1) 

W2(x;a,6,c) = —- [ y^"^(y - 1)"“(2/- dy 

r(a - c +l)r(l - a) 

{x G (—OO, 0) U (0, 1), Re (c — 1) < Re a < 1), (7-2) 

W3ix;a,b,c) = f \y\^~^\y - M~ly - dy 

T{a- c + l)T[c-b) Jy/x>i 

{x G (—OO, 0) U (1, oo). Re (a — c + 1) >0, Re {c — b) > 0), (7-3) 

wYx;a,b,c) = y'’"^(l - y)"“|x - dy 

{x G (—oo, 0) U (1, oo). Re a < 0, Re 5 > — 1), (7.4) 

W 5 (x; a, b, c) = - y^ix - yT~''~^ dy 

r(a-c + l)r(5) 

{x G (0, oo). Re 6 > 0, Re (a — c + 1) > 0), (7.5) 

m6(x;a,6,c)= ^|"~^ ~^^^ x^-^ [ yb-^l - y)-\x - yr^-^y 

r(l-a)r(c-6) Jo<^<i 

{x G (0,1) U (1, oo). Re (c — 6) > 0, Re a < 1). 


(7.6) 
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Euler type integral transformations with integrand \yY 


y 


\c—b—l 


X 


y\ 


u;i(x;a, 6 ,c) = '^{y - ^ ^{y - x)' 


W 2 {x-,a,b,c) = 

W 3 {x;a,b,c) = 

Wi{x; a, b, c) = 

W 5 {x-,a,b,c) = 

wq{x; a, b, c) = 


T{b)T{c-b) 

{x < 1, Rec > Re 6 > 0), 
r(2 - c) f 


dy 


\y\--(l-yy-b-^\x-y\-'^dy 


T{a-c + l)r(l - a) 

(x € (—oo, 0) U (0,1), 2 > Re (a + 1) > Re c), 
r(a -6 + l) 


T{a-c + l)T{c-b) Jo ^ ^ \ y\ y 

{x G (—oo, 0) U (1, oo), Re (a + 1) > Re c > Re 6), 

r(6-a + 1 ) 


[ \y\--\y-l\c-b-l\y_^\-ady 

J v/x>l 


r(—o)r(6 + 1 ) jy/x>i 

{x G (—oo, 0) U (1, oo). Re a < 1, Re 6 > 0), 
r(a + 6 -c + l) 


[ {-yy--{l-yr-^-\x-y)-^dy 

J —oo 


r(a-c + i)r( 6 ) j_, 

{x G (0, oo), Re 6 > 0, Re (a + 1) > Rec), 
r(c — a — 6 + 1 ) 


r(i-a)r(c-6)yo<^<i ' ' ' I ^ 

{x G (0,1) U (1, oo), Re (c — 6 ) > 0, Re a < 1). 


(7.7) 

(7.8) 

(7.9) 

(7.10) 

(7.11) 

(7.12) 


The integrals in (7.1), (7.3), (7.6) (for tri, W 3 , wq) and in (7.8), (7.10), (7.12) (for W 2 , W 4 , wq) 
are of fractional integral type. Formulas (7.3), (7.8), (7.10) and (7.12) can be reduced to (1.1) 
by a change of integration variable, while (7.1) is a rewritten version of (1.1). Each of the six 
formulas (7.1)-(7.6) is paired with one of the six formulas (7.7)-(7.12) by a transformation of 
integration variable which involves x. The pairing is as follows: 


(7.1) GG (7.9), (7.2) GG (7.10), (7.3) gg (7.7), 

(7.4) GG (7.8), (7.5) GG (7.11), (7.6) gg (7.12). 


In 1874 Letnikov [14, (11), (13), (16), (17), (21), (22)] (see also Sostak [23, Section 4]) already 
gave the Euler type integral representations (7.8), (7.12), (7.10), (7.1), (7.6), (7.3), respectively, 
as solutions of the hypergeometric differential equation [14, p. 115, (A)]. In order to arrive 
at these results he essentially considered the fractional integral operator as a transmutation 
operator with respect to the hypergeometric differential operator. 


int. rep. 

transform 

from 

to 

case 

(7.1) 

(5.1) 

rci 

rci 

c+ 

(7.2) 

(6.4) 

W3 

W 2 

o+, c+ 

(7.3) 

(5.3) 

W3 

W3 

c+ 

(7.4) 

(6.3) 

Wl 

W 4 

fl+, c+ 

(7.5) 

( 6 . 1 ) 

Wl 

We 

c+ 

(7.6) 

(5.6) 

We 

We 

c+ 

(7.7) 

( 6 . 20 ) 

W 4 

Wl 

a—, c— 

(7.8) 

(5.7) 

W 2 

W 2 

a— 

(7.9) 

( 6 . 10 ) 

We 

W3 

a— 

(7.10) 

(5.8) 

W 4 

W 4 

a—, 6 —, c— 

(7.11) 

(6.9) 

W 4 

We 

a— 

(7.12) 

(5.9) 

We 

We 

a—, b—, c— 
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The above table gives for each Euler type integral representation in the first column the 
transformation formula in the second column from which it can be obtained by specialization 
of parameters (6 = c in the first six rows and a = — 1 in the last six rows). The transformation 
formula sends Wi in the third column to wj in the fourth column. Its case is given in the fifth 
column. Formula (7.5) is essentially the same as [8, 14.4(9)]. 


8 Generalized Stieltjes transform 

and fractional integral transform combined 


As observed in [8, p. 213], generalized Stieltjes transforms of different order are connected with 
each other by fractional integration. Formula (8.1) below was essentially given there, and later 
given with proof in [10, Theorem 9]. Some related identities can also be proved: 


Proposition 8.1. Let Re(l — z^) > Re/i > 0. Assume that f G L\^^{{m,M)) and that the 
integrals on the right-hand side of the four identities below converge absolutely. Then 


r if 


f 


rM 


\z-yy- 

r(i-^) 


rM 


fiz) 


’ m 

ry 


{y - z) 


l-u 


Ax-yY ^ 

' r(rt = 

dz\iK^dy = 




T{l-u) 


(x<m),(8.1) 

(M<x),(8.2) 


dy = j^ (x<m), (8.3) 

. - dy = [ f{z)-r^ - ^dz (x > M).(8.4) 

(2/-x)i- J-oc A-zy-^^- ^ ’ 


{x-yY '' 
T{l-u) 


The proofs are immediate, by the Fubini theorem and by a version [18, (5.12.3)] of the beta 
integral. Furthermore, (8.2) is an immediate consequence of (8.1), and similarly (8.4) of (8.3). 

Examples of formulas (8.1)-(8.4) can be found by combining suitable fractional integral for¬ 
mulas in Section 5 with suitable generalized Stieltjes transform formulas in Section 6. For 
instance: 


• In (8.1) let m = 1, M = oo, f{z) = z'^~^WG{z;a,b,c) and use (5.2) and (6.2). 

• In (8.2) let m = —oo, M = 0, f{z) = {—zY~^wi{z;a,b,c) and use (5.5) and (6.1). 

• In (8.3) let m = 1, f{z) = z‘^~^we,{z;a,b,c) and use (5.6) and (6.2). 

• In (8.4) let M = 0, f{z) = {—zY~^wi{z;a,b,c) and use (5.1) and (6.1). 
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